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The lensing effect of a cosmic string is studied, and some new methods are proposed to detect 
the cosmic string. The technique for using jets as extended gravitational lensing probes was 
firstly explored by Kronberg. We use the "alignment-breaking parameter" rjc as a sensitive 
indicator of gravitational distortion by a wiggly cosmic string. Then, we applied the non- 
constant deflection angle to jets, and ijq of a specific jet is just related to the projected slope of 
the jet. At least three jets in the sample of Square Kilometer Array (SKA) would have significant 
signals (rjc > 10°) if the wiggly infinite cosmic string existed. The distortion of elliptical object 
is also studied and used to do a statistical research on directions of axes and ellipticities of 
galaxies. In the direction of the string, we find that galaxies appear to be more elliptical for 
an observer and the distribution of apparent ellipticity changes correspondingly. Ellipticity 
distribution of current SDSS spiral sample has the signal-to-noise ratio up to 8.48 which is 
large enough for astronomical observations. The future survey, such as Large Synoptic Survey 
Telescope (LSST) and Dark Energy Survey (DES), would weaken the requirement of special 
geometry in the data processing. As a result, all kinds of distributions, including ellipticity 
axis distribution, would serve as probes to detect wiggly strings in the near future. In brief, 
if a wiggly cosmic string existed, these signals would be convenient to be observed with the 
future weak lensing survey or other surveys in deep space. If there was no lensing signal in 
these distributions, it would give the upper limit of the abundance of infinite strings. 

Keywords gravitational lensing, jets and outflows, cosmic string, galaxy 

PACS numbers 98.80.Cq, 98.62.-g, 98.62.Nx, 95.30.Sf 



1 Introduction 

A cosmic string, as a kind of topological defect, is gen- 
erated by symmetry breaking of unified theories, such as 
fundamental string theory, brane inflation, M-theory, etc 
[U [21 13] . The detection of a cosmic string can be used 
to constrain different inflation models because strings 
were formed in the symmetry-breaking phase transition 
of some special inflation models, such as brane inflation 
[31 E] . Cosmic string gas may produce observable non- 
Gaussianity if the string scale is at TeV scale [7]. Ad- 
ditionally, as a special kind of "dark matter", a cosmic 
string can produce cosmic ray positrons by the cosmic 
string cusp annihilation process, which might explain the 
positron excess of PAMELA and ATIC experiment [8]. 
Moreover, a cosmic string may play a sub-dominant role 



in the formation of large structure of the universe to some 
extent OHO]. 

The gravitational properties of strings are studied with 
the weak-field approximation. According to the metric 
of a string, the light is deflected by an angle of AnGU . 
Thus, this gravitational effect could give rise to double 
images of background cosmic objects within the angle 
of order 5<p from the string [TU [121 HHj- Their gravi- 
tational effect is measured by a dimensionless parameter 
GU/c 2 , which is about 1CT 6 for GUT (grand unified the- 
ory) strings. Due to the gravitational effect caused by a 
cosmic string, many observations can be used to probe 
the properties of a cosmic string. The current result is 
GU/c 2 < 0.7 x 1(T 6 , which is obtained by CMB con- 
straints [13] • First LIGO search for gravitational wave 
bursts from cosmic strings also give nearly the same up- 
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per limit [15]. A cosmic string makes the space-time 
around it to be a conical space-time, resulting in the 
appearance of undistorted double images [T|)J H7J ITS] , 
Micro-lensing provided a way to detect lensing signal 
even when the image splitting is too small to resolve with 
astronomical measurements [19] . With the development 
of the technology of weak lensing, the weak lensing effect 
produced by a cosmic string could also be a useful tool 
to constrain the properties of a cosmic string |20L I21j . 

In this paper, we consider the weak lensing effect 
caused by a cosmic string, and polarized jets and ellipti- 
cal objects are used as the background sources. Using the 
fact that the polarization angle of a background radio jet 
is not changed by gravitational distortion [22 , it is con- 
venient to define an " alignment-breaking parameter" t/q 
[23], which is sensitive to gravitational distortion. This 
method was used to analyze the polarization of the jet 
3C 9 and give the redshift and properties of lens galax- 
ies [23]. Theoretical modeling of weakly lensed polar- 
ized sources is developed [25] for the prediction of effects 
of gravitational lensing. As promising as it is, gravita- 
tional lensing is hampered by the lack of the knowledge 
of the intrinsic morphology of the sources being lensed. 
However, the polarization of the source could supply a 
useful probe to determine the intrinsic morphology of 
the source. Nevertheless, the lens employed by them are 
all spherical galaxies or elliptical ones. To our knowl- 
edge, this is the first investigation of weak lensing ef- 
fect of a cosmic string with polarized sources as back- 
ground objects. Once the background galaxy was lensed 
by a wiggly cosmic string, the ellipticity axis distribu- 
tion (EAD), and the ellipticity distribution (ED) would 
be changed. Though there might be just several infinite 
straight strings [18] . they would distort the light pass- 
ing by it with any impact parameter. For a straight 
string, the deflection angle of an image is constant. Ac- 
tually, there is no absolutely straight string. And the 
wiggly string would produce non-constant deflection an- 
gle which would result in the distortion of background 
sources [20 . In this paper, we will consider different kinds 
of sources as background objects to probe the property 
of the wiggly cosmic string. 

The paper is organized as follows: First, three ba- 
sic concepts, alignment-breaking parameter, EAD and 
ED, are introduced and the deflection angle is deduced. 
Second, we use the "alignment-breaking-angle" to study 
lensing effect of a straight jet. EAD and ED of galaxies 
are also studied. Then, we put the formulas into sim- 
ulations and provide the strategy to probe the signals 
produced by weak lensing effects of a cosmic string. In 
the part of conclusion, we propose some statistical meth- 
ods to search for such signals in the future surveys. 

2 Theory 



2.1 Basic concepts 

2.1.1 Alignment-breaking parameter rjc 

In order to measure the signals yielded by a cosmic 
string, we should at first define a parameter rjc |23| as 
the angle between the direction of observed polarization 
and the tangential vector of a jet: 

VG (6) = yj(ff) - xo(ff) + <0) (1) 

where ifj(9) is the angle of tangent vector relative to the 
reference line of the jet at the projected position 6. xo is 
the polarization vector, and n is some intrinsic deviation 
from a perfect jet (see Figjl]). Define e img ,e p , e real ,e ref 
as the unit tangent vector of the fiducial line of the im- 
age, direction of polarization, fiducial line of the source, 
direction of reference line, and the vectors corresponding 
to ip,Xo, k ar 

Xo =e p — e re f, 

VG &i?ng ^reah (^) 

so rje = tp + n — Xo- Vg is a measure of the projection of 
shear onto the jet itself. 

In equation ([2]), the polarization, xoj tangent angle, 
ip, and the intrinsic deviation of the polarization from 
the tangent angle, kq, must be measured along the jet or 
lobe. However, the intrinsic polarization of the source, 
Xo, is coupled with the Faraday Rotation, which causes 
the polarization of a photon to rotate as it traverses a 
magneto-ionic medium. By measuring Faraday rotation 
at different wave lengths, one could eliminate the varia- 
tion of rjc brought by Faraday rotation. The true value 
of rotation measure is needed if we want to rule out the 
constant part produced by Faraday rotation [26) . In ad- 
dition, the value of k is a source-dependent angle that 
allows for a source-intrinsic variation of t\q. Fortunately, 
the polarization of a jet is highly aligned with the local 
direction of the jet [371 HH HH] • Simulations of the intrin- 
sic magnetic field of a radio jet are conducted by using 
MHD model of a jet [27J. Constrains on k would re- 
quire the observation of a low-redshift sample of sources, 
and even then, k could only be understood in a statis- 
tical sense [26] , Therefore, with the mature theory of 
the magnetic field of jets, we could independently derive 
the value of k(9) and consequently use t\g to reveal the 
strength of the gravitational bending varies along the jet. 
That is, observing a jet and its polarization could confirm 
the value of rjQ, which is deduced from the gravitational 
effect of a cosmic string. Thus, t\g serves as another 
parameter to reveal the property of a cosmic string. 



Fa-Bo Feng, Radio Jets and Galaxies as Cosmic String Probes, Front. Phys. 



111 




Fig. 1 The relation of the alignment-breaking parameter, t\q, to other angles k, ip, xo- The dotted (dashed) line 
represents polarization (image) , which is the projection of the direction of polarization on the image plane . 



2.1.2 The ellipticity and ellipticity axis angle of the 
galaxy 

For a spiral or elliptical galaxy, the intrinsic shape might 
be simplified into a triaxial ellipsoid, x 2 /a 2 + y 2 /b 2 + 
z 2 /c 2 — 1 with a ^ b ^ c. The galaxies can be divided 
into three kinds: oblate galaxy (a=b), prolate galaxy 
(b=c), and triaxial galaxy (a =/= b =/= c). Define 7 = c/a 
as the thickness of a galaxy and e = 1 — b/a as the in- 
trinsic ellipticity, we can convert the intrinsic coordinate 
system to the observer's coordinate system and compute 
the resulting apparent axis ratio e = b b/a t, and the 
ellipticity axis angle (EAA) ^ of a triaxial galaxy [30] : 



X + Z-^{X-Z) 2 +Y 2 
X + Z + ^/(X-Z) 2 +Y 2 



1/2 



VP = - arctan 
2 



Y 



X - Z 



(3) 
(4) 



where 



X = [1 - e(2 - e) sin 2 ip] cos 2 9 + j 2 sin 2 
Y = e(2 — e) cos 9 sin 2ip, 
Z = l-e(2-e) cos 2 tp, 



and 9, cj> is the spherical coordinates of the line of sight. 
The photometric analysis of Sloan Digital Sky Survey 
(SDSS) provide a model-free measures of the axis ratio 
by the technique of adaptive moments. This method 
guarantees the accuracy of the values of a b,b b, which 
are the long and short semi-major axis lengths. By fit- 
ting models to the shape measurement of SDSS Data Re- 
lease 1 (DR1), the observed apparent axis ratios could 
be modeled by adopting a Gaussian distribution of 7 and 



a log-normal distribution of e [31] : 



/(7) oc exp 



(7-M7) 2 



2a 2 



1/2 



/(e) cx - exp 



(lne-^) 2 
2a 2 



< 7 < 1, (5) 
lne<0. (6) 



With four parameters /i 7 , cr 7 , /i, a estimated from the 
data, we could describe the shapes of elliptical galax- 
ies and spiral galaxies with the distribution of apparent 
axis ratio. 



2.2 Deflection angle by a wiggly cosmic string 

We quote the result of ref. [32] for the deflection of a light 
ray by a static wiggly cosmic string, 

rD s -D L 

a = -V_L/ $(x l ,x 2 ,x 3 )dx 3 , (7) 

J-Dt 

where Dg and Dt, are separately the distance from the 
observer to the source and to the cosmic string. Notice 
that the analytic expression of the deflection angle does 
not include the constant deflection angle which is of no 
use for distorted background objects. The deflecting po- 
tential is given by 



$>(xi,x 2 ,x 3 ) = 2G 



(8) 



Here G is the gravitational constant, x is a reference 
point, = dx^/dX is the tangent vector to a geodesic 
g which is chosen as reference (see Fig{2]). Also, r(£) is 
a parametrization of the string, in terms of which 



b 



(9) 
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where T^ v is the stress-energy tensor and r/^ is the 

Minkowski metric, i.e. rj^ = diag( h H — h). Generally 

speaking, the stress-energy tensor of a cosmic string is of 
the form 

Ti iv = Vu n u v _ tw, ( 10 ) 

where U is the energy per unit length, T is the tension 
which doesn't equal to U for a wiggly string, it M is the 
4-velocity of the observer (w^ is also a timelike vector, 
i.e. u^u^ = — 1) and is a spacelike {v^v v = 1) unit 
vector tangent to the string word sheet. Let us consider 
a tilted static cosmic string aligned as Fig(2] illustrated, 
and the two unit vector are 



Then, we recall the form of deflection angle in ^ , which 
reduces to 

al =a2 

4G(T - U) {2x 2 A0A4> + Xl (AO 2 - A</> 2 )) D s 
~ {AO 2 + A4> 2 ) 2 D L (D L -D S ) 

(15) 

From the equations above, we find that the deflection 
angle is proportional to {AO 2 + A^ 2 )^ 1 approximately. 
To compare the deflection angle deduced from ([8| with 
that deduced from ( 13 ), we have plotted the graphs of a 



^=(1,0,0,0), 

=(0, — sin sine/), cos <f>, — cos sin < 



(11) 



We take the tangent vector of the reference geodesic to 
be W = (1,0, 1,0). From we deduce that 



,9 =2^}? 

=U -Tsin 2 (bcos 2 0. 



after adopting the above approximation(see Figj3|. As 
a result, we find that the general shapes of them are 
similar with those of Figj3] However, the sign of the 
two deflection angle is not opposite now. Because we 
take Xi/Dg — x 2 /Dg — 10~ 6 and the deflection angle 
is up to 10 -7 , the distortions of the background objects 
are large enough to be observed. The statistics of back- 
ground jets and galaxies would reveal the secrets of the 
wiggly cosmic string. 



(12) 3 Simulations 



If we take = 0, then & = U — T sin 2 0, which is the re- 
sult of ref. [32 . When U{£) and T(l) are constants, the 
deflecting potential, after integration over X3, reduces to 



x 2 , x 3 ) = - AG{U - T sin 2 0cos 2 0) ln(f 2 - 

(— x\ sin sin <f> + x 2 cos <j> — x% cos 9 sin 0) z 

(13) 

In the limit where \x$\ ^> | cci | , j ara ] and under the con- 
dition that {0,(f>} is not very close to {0,7r/2}, the de- 
flection angle, i.e. a(xi,X2) = {al,a2}, is deduced from 
However, the form of it is complicated, so we have 
plotted the 3D graph of a{0, </>) (see Fig|3| after taking 
GU = 1.1 x 10" 6 , GT = 0.9 x 10" 6 . Jleasonably, we take 
p = D L /D S = 1/2 and , \xi\/D s = \x 2 \/D s = lO" 6 . 

Illustrated by Fig(3j the deflect ion angles, al and a2, 
are sensitive to the change of {0, </>} especially when they 
come close to {0, 7r/2}. Additionally, the signs of the two 
deflection angles are almost opposite. However, this re- 
sult is limited by its assumption of the value range of 
{0, 4>}. Therefore, it is necessary to explore the form 
of deflection angle in the limit where {0,(j>} — > {0,7r/2} 
in order to enhance the value of deflection angle. After 
definition of {0, cf)} = {A0,tt/2 + Acf)} and the limita- 
tion where AO — > and A(j> — > 0, the form of defecting 
potential of ^ reduces to 

$(x u x 2 ,x 3 ) = 4G{U-T)\n[(x 3 A0+x 2 ) 2 +(x 3 AO+ Xl ) 2 }, 

(14) 

where AO and A(f> should be larger than 10 6 which 
makes x 3 A0 and x 3 A(f> much larger than x\ and x 2 . 



3.1 Radio jets 

Considering the statistical errors and systematic errors, 
we can only extract the signal of r]c > 10° (i.e., t\g > 
0.175 rad) |26j from the polarization analyzing of ra- 
dio jets. So we should take the approximation that 
'AO, A(j> — > to yield large distortion angle. 

First, we have parameterized the configuration of the 
radio jet as follows: 



xi =r cos ft + x\ c , 
x 2 =r sin f3 + x 2c , 



(16) 



where r is the distance from reference point to the cen- 
ter of the jet, and (3 is the angle between its vector 
and the axis of X\. The distortion could be deduced as 
x[ = xi+Dsal and x' 2 = x 2 + D$a2. The corresponding 
r\G reduces to 

VG =AP 

= arctan[(4G(T - U){A0 2 - Acf 2 ) + (8G(T - U)A0Acj) 
+ P (A0 2 + A4> 2 ) 2 - P 2 {A0 2 + A4> 2 ) 2 ) tan/3) 
/(4G(T - U){A0 2 - A(j) 2 ) - (-1 +p)p{A0 2 + Acf 2 ) 2 
+ 8G(T - U)A0Acj) tan /?)] - (3. 

(17) 

From the equation above, we find that r\Q is indepen- 
dent with the coordinates of the points on the jet, i.e. 
r\Q is constant for a specific straight jet. However, r\ G 
is sensitive to the inclined angle /3. The shape of t\q is 
different with the sinusoidal r\Q shape yielded by spher- 
ical lens [23]. However, cosmological birefringence could 



Here, G, U, T are all in nature unit. We adopt the value of GUT string which is produced during the GUT transition, E = lC^SGeV . 
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Fig. 2 The wiggly string leads to some distortion of the straight jet. 




Fig. 3 The left graph shows the 3D plot of al and the middle one shows that of a2. Notice that the two graphs 
are all plotted in the limitation where {9,cf>} is not close to {0, 7r/2}. The right one illustrate the value of a when 
{61, (/>} -> {0, tt/2}, where {0, 0} = {A0, vr/2 + A0}.. 



1o 

1" I 




Fig. 4 The left graph shows the 2D plot of tjg and the right one shows the log plot of rjc as a function of ft. Notice 
that the two graphs are all plotted by taking integer k from to 40 colored from red to blue. 
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Fig. 5 The distribution of t\g is illustrated by this graph in which log-log coordinates make the graph more 
appropriate. There are about 100,000 points for statistics. 



also result in a rotation of the polarization of photons, 
which could yield a nonzero 7]g [33 . Though this kind 
of effect may appear similar with that caused by a cos- 
mic string, they have different features. Birefringence is 
a cosmological effect that could be detected in the whole 
universe with nearly the same order of magnitude. That 
is, if t]g 7^ for one jet, it would also be founded in 
other jets. However, the signature given by a wiggly 
cosmic string is local. Of course, other sources, such as 
gravitational waves and rotational black holes may also 
yield nonzero rjQ, but their shapes are not constant for 
different points on the jet [HI [M]. Additionally, with cur- 
rent resolution of the radio telescope, such as SKA, the 
main errors come from the intrinsic noise of the polar- 
ization k, the error of choosing the anchor points, and 
the residual noise after extracting the effect of Faraday 
Rotation measurement (RM). Considering these factors, 
we can only extract the signal of t\g > 10° [33] from the 
polarization analyzing of radio jets. Thus, it is impor- 
tant to yield large enough signals by carefully tune the 
parameters. With current observation, the properties of 
the cosmic string are parameterized by GU = 1.1 x 10 -6 
and GT = 0.9 x 10~ 6 ; and the Dl to D$ ratio p is taken 
to be 1/2. We could plot the 2D graph series with differ- 
ent ranges of AO and A(j) (see Fig. [4j. From these two 
graphs, we find that there would be signals large enough 
to be observed when AO = A<p ~ 0.001. Notice that the 
value of the signal is always proportionate to G(U — T). 
And rjQ is sensitive to the angle of inclination /3, which 
also give us the information of the configuration of the 
jet. In order to get a general idea of the probability for 
the detection of a cosmic string through analyzing t]g, we 



put the distribution of tjg into programming and draw 
the graph by Mathematica. AO, Aa are randomly chosen 
in the range of [—0.005,0.005]. Larger interval result in 
a little proportion of useful points though the sample is 
large at the interval; and a smaller interval would lead to 
a very small sample which is not appropriate for statis- 
tics. The number of background objects was taken to be 
100,000, which calls for 100,000 times of random choices 
of these parameters in their intervals. The parameters of 
cosmic string are set to be the same as the above ones, 
and the lower limit of rotation angle or t\g is taken to 
be t\g > 10° (i.e., r\G > 0.175 rad). After that, the final 
result is revealed by Fig(5j 

In this histogram, we find about 1/5 jets would yield 
rotation angles larger than 10° when AO and Aa are at 
the interval of [-0.005,0.005]. Notice that 100,000 jets 
in such a small interval do not accord with reality. We 
adopt such a large sample is just for getting the proba- 
bility of producing large signals. Though, the population 
in such an interval is small currently, SKA may change 
the situation. It would carry out an All-Sky RM Survey, 
which could yield RMs for approximately 2 x 10 4 pulsars 
and 2 x 10 7 compact polarized extragalactic sources in 
about a year of observing time [35] . In this All-Sky RM 
survey, we may get useful information of the intrinsic po- 
larization angle of pulsars who have jets and lobs. To be 
on the safe side, we consider 2 x 10 6 radio objects as an 
ideal sample. If the wiggly cosmic string or the string 
web exists, the population of the useful jets in the sub- 
sample of the right area is 1/5x(0.01) 2 x2x10 6 /(4tt) ~ 3. 
And in other areas of the sky where AO and Acf> larger 
than 0.01, there would be nearly no observably signal. 
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Notice that these jets are characterized by constant ro- 
tational angle t\q. 

In sum, the value of t\g is determined by the relative 
position of the jet with respective to the cosmic string 
and the inclined angle of the string. The position where 
r\Q would be enhanced is in the direction of the string; 
and the larger the inclined angle is, the larger the deflec- 
tion angle is. Additionally, constant r\Q resulted from the 
straight jet is easily to be extracted from other noises. 
If t\q is irregular, it would give us information about the 
configuration of the jet. In the near future, with the 
improvement of the accuracy of RM measurement and 
the maturation of theoretical model of the intrinsic mag- 
netic field in the jet, r\Q would show its strong power of 
detecting the wiggly string. 

3.2 Galaxies 

In order to study the cosmic string with current large 
scale survey, it is best to use galaxies as background 
sources. And the weak gravitational effect from a cosmic 
string should be analyzed with statistical method. Ad- 
ditionally, the method used here is a little different from 
that used in the statistics for the cosmic shear. The fol- 
lowing method only abstracts the local information, and 
the study of cosmic shear aims at getting statistical infor- 
mation. The ellipticity axis distribution (EAD) of circu- 
lar 2D sources was used to examine the shear signature 
of nearly circular sources by Dyda and Brandenberger 
[2"U] , In our paper, EAD might also be a good parame- 
ter to describe the distortion of the EAA of an arbitrary 
galaxy in the data processing of the future survey. Ad- 
ditionally, we used the ellipticity distribution (ED) as a 
probe for the cosmic string detection. Rather than us- 
ing the circular object as the source, we extended the 
shape of the source to be a 3D triaxial object. Then we 
project the 3D image onto the surface perpendicular to 
the viewer. Points on the surface of the elliptical object 
have coordinates: [x\ c -I- xi, x^ c + x%, x^ c + X3). 

In the following formulas, we use e = b b/a b to be the 
apparent axis ratio. And the luminosity density I is 



I = l(a v ), 



ihere 




(18) 



The above formula is used by Stark and Binney to 
calculate the apparent axis ratio [351 [30]. However, 
we adopted a different transformation from the two- 
parameter transformation by them. They just use 9, <j) 
to fix the system of the observer by making x' axis lie 
in the O-xy plane j30j. we adopted a, /3, 7 to describe 
the system of the observer. The system of {xq, yo, z Q } is 
a system in which the axes of Xq , yo , z are aligned with 
the principal body axes of the ellipsoid. After transfor- 
mation from this system to the coordinate system of the 



observer, their coordinates would reduce to 
xq =x\ cos P cos 7 — X3 sin /3 + X2 cos j5 sin 7, 
Ho — x 3 cos P sin a + X\ (cos 7 sin a sin j3 — cos a sin 7) 

+x 2 (cos a cos 7 + sin a sin /3 sin 7) , 
z o = x 3 cos a cos P + x\ (cos a cos 7 sin /3 + sin a sin 7) 

+X2 (— cos 7 sin a + cos a sin /3 sin 7) . 

With the above transformation, we could get the el- 
liptical radius, a v [x\, X2, X3), in the coordinate system 
of observer. With the method explored by Stark and 
Binney [30J, the apparent axis ratio of the ellip- 
soid on the sky will be easily deduced. Due to the 
complexity of formulas, the calculations are left for the 
computer to do. At last, the change of the ellipticity 
distribution (ED) is deduced from the distortion of e, 
Ae = e(a, (3, 7, C, 9, </>) — e(a, (3, 7, C, £)■ Therefore, both 
ED and EAD are determined by distributions of these 
parameters. Among these parameters, the distributions 
of a,/3,7 are determined by the random distribution of 
points on an unit sphere; 9 and <p are replaced by A9,A<f>, 
which obey the flat distribution; and (, £ obey Normal 
and Log-Normal distributions respectively. 

It is obvious that Ae and t\q are determined by the 
value of a. That is, inserting (15) in a v [x\ + Dsa,\,X2 + 



DsU2,xz), we deduce the distorted apparent axis ratio 
(|3| and EAA Q. Ae would also contribute to the cos- 
mic shear if a wiggly cosmic string existed. The shear 
signal yielded by the cosmic string is subtle when 9 -/> 
and <j) tt/2. However, it may contribute a little frac- 
tion to the cosmic shear and result in systematic errors. 
The shear field of the wiggly cosmic string should also be 
studied for this reason [3TK37]. If A6, Acj) -> 0, Ae would 
be proportional to (A6> 2 + A</> 2 ) -1 and increase even to 
the order of 0.1. It will be interesting to find these sta- 
tistical signals produced by the network of strings in the 
future weak lensing survey. 

Similarly, with the statistical method used in the last 
section, we could get the variation of signal-to-noise ra- 
tio (SNR) with the length of the interval of A9 and A<p, 
i.e. [-6/2,6/2] = [A9 mm ,A9 max ] = [-O.OOlfc, O.OOlfc]. 
It is necessary to study SNRs of different modes (spec- 
ified by integer k) of sky survey to know the feasibility 
of this method. Because the signal of ED/EAD and the 
noise of the sample in the area of S 2 increase with 5~ 2 
simultaneously, we should find the optimal 6, i.e. 6 opt , to 
make SNR highest. The specific processes are as follows: 

First, considering ED/EAD without a string (origi- 
nal ED/EAD), we produced a realization of ellipsoids 
according to the distributions explored by Padilla and 
Strauss [38]. Here, we assume the empirical ED is 
ideal because the fitting error is less than the statisti- 
cal error especially in the interval of apparent axis ra- 
tio, e £ [0,0.1]. The population of the sample shall 
be calculated according to the surface number density 



Vlll 
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of the survey sample. For example, there are about 90 
targets per square degree in the main galaxy sample of 
SDSS DR6. Thus, a realization of ED/EAD would call 
for about 90 S 2 x 180 2 /7r 2 times of generations of ran- 
dom points. The parameters included in ED/EAD are 
a, 0, 7, £, £. a,/3,7 are determined by the random dis- 
tributed points on an unit sphere. 7, £ are chosen ac- 
cording to the empirical ED. A9 and A</> are randomly 
chosen in the interval of [—O.OOlfc, O.OOlfc]. 

Second, by repeating the first step 100 times, we com- 
puted the mean and standard deviation of the expected 
number of original ellipsoids in each of the 10 bins. Ad- 
ditionally, the mean of distorted ellipsoids is also gotten 
simultaneously. The standard deviation is the function 
of S. By the way, the reason for setting 10 bins is that 
the population of every bin is limited by 0.002k. Too 
many bins would result in large noise in every bin. At 
the same time, through this process, we have produced a 
large sample of ellipsoids distorted by the string because 
it is 100 times more than a single step. Such a large 
sample would realize ED and EAD of the lensed galaxies 
with small statistical errors. 

Third, subtracting the original ED/EAD from the 
lensed ED/EAD, we extracted the useful signal. Then, 
divide the signal by the standard deviation of the original 
ED/EAD, we got SNR as a function of S (Fig|6]). 

We adopted two different ED models, one of spirals 
and the other of ellipticals, and computed SNRs of them 
respectively. Due to the little contribution of fitting er- 
rors and systematic errors to the total error in the first 
bin, only statistical errors would pollute the signal. So 
we divided the signal by the statistical error to get SNR. 
After computing 34 SNR graphs (17 for EDs and 17 for 
EADs), we found that SNR of ED in the first bin, e < 0.1, 
is much larger than that of the others. However, in most 
of the cases, the SNR of EAD is less than 1. Thus EAD 
could not be used to detect the cosmic string as Dyda 
and Brandenberger predicted in their paper [20] . Thus, 
the SNR of ED in the first bin is chosen to mark the 
level of feasibility. In addition, we chose two surface 
number densities, one of which is the density of SDSS 
main galaxy sample and one is the lower limit of the 
future LSST's density. The types of galaxies also con- 
tribute to the variation of EDs and EADs. Thus, we have 
four optimal modes corresponding to four distributions: 
2 types of galaxiesx2 kinds of number densities. Each 
mode require one specific survey. For example, it would 
be better to divide the sky area into small circles, each of 
which has an area of S 2 . Notice that every circle is sup- 
posed to have coinciding area; and every inscribed square 
of each circle should adjacent to those of surrounding cir- 
cles. However, if the sample is large enough, the require- 
ment of specific geometry might be weak. We could find, 
in Figjm four optimal modes, corresponding to 5 op t) for 



the analysis of the current and future samples. 

In Fig|6) the re are four SNRs for four EDs, 
each of which is corresponding to one specific den- 
sity and one type of galaxy. It is easy to 
get the optimal modes of each SNR. The optimal 
modes for 90 spirals / degree 2 , 1000 spirals / degree 2 , 
90 ellipticals /degree 2 , 1000 ellipticals /degree 2 are 
A9 opt = 0.01, 0.004, 0.008, 0.002. The optimal SNR 
of spirals with current number density is 8.48 which is 
larger than that of ellipticals 3.35. With higher density, 
ellipticals would also be appropriate for abstracting lens- 
ing signals because of its high SNR 18.53. The results 
indicate that the mode of spirals is better to be used as 
a probe, particularly for current samples with low num- 
ber density of 90 spirals /degree 2 . Though spirals suffer 
more extinction of dust than ellipticals do, the ED of 
them is still the best probe for cosmic string for its high 
SNR. Besides, the development of the correction for dust 
extinction and the division of the sample according to 
magnitudes would improve the accuracy of this method 
[38] . For the future samples with higher number density, 
it is better to adopt all available modes (with different 
S) rather than only one mode. By the way, the large 
amount of data would decrease the statistic errors. Un- 
der this condition, the systematic error might be larger 
than the statistic error. The development of this method 
would require accurate models of EDs for different types 
of galaxies to decrease systematic errors. 

However, the SNR of EAD is too low for usage. Thus 
we only give the error bar of the EAD with respect to 
each optimal mode. There is almost no SNR larger than 
1. The result indicates that EAD has relatively low SNR 
and is not so appropriate for detecting cosmic string with 
current data. The low number density of the sample im- 
pedes the development of this method. However, EAD is 
independent with ED models because the original EAD is 
always flat. Besides, the release of the future data would 
make the statistic errors lower than systematic errors. 
A density of 10,000 galaxies /degree 2 might be enough 
for EAD to be used to abstract useful lensing signals. 
Thus, EAD would even be better than ED in a long run. 
In addition, for the reason that the future large sam- 
ple shall weaken the requirement of optimal modes, the 
usage of different modes simultaneously would of course 
increase the power of these methods. Nevertheless, what- 
ever modes we use, ED and EAD would be abnormal if 
the galaxies lied in the direction of the string. The po- 
sition of these abnormal EDs and the special angle of 
EADs would reveal the true direction of the string. The 
reason is that the change of the population in the middle 
bin of EAD, i.e. \& = 0, is larger than any other bins. 
This special direction is parallel with xl axis, which in- 
dicate the direction of the cosmic string. 

Generally, these methods are different from the one 
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developed by Dyda. Dyda and Brandenberger employed 
special range of the ellipticity to get obvious EAD [2D] . 
Here, we have used special position to get abnormal EAD 
and ED. It is easier to find specific area than to find spe- 
cific ellipticity because the improvement of survey area is 
easier than the correction of ellipticity. Moreover, the fu- 
ture survey will weaken the requirement of survey modes 
for its lower SNR. For example, the optimal ED mode of 
ellipticals with high number density is 5 op t — 0.002. Un- 
der this condition, we can just get a small solid angle of 
(5 2 180 2 /7r 2 « 0.013 deg 2 . The survey area of LSST will 
include 30,000 deg 2 in which a databases of 10 billion 
galaxies would be produced [39] ■ Thus, we can easily get 
4000 galaxies in 0.013 deg 2 which would largely decrease 
the statistic noise and increase SNR greatly. However, 
if there were no significant variations of ED and EAD, 
we would put stronger constraints on string qualities, 
such as G(U — T). In (14), a is in inverse proportion 
to G(U — T). That is, there would be two factors con- 
tributing to the result that no ED signals were detected. 
The first one is G(U~T) < 10~ 6 , and the second one is 
that the area covered by the survey is not large enough. 
With larger and deeper sky areas being surveyed, the 
first factor would be the main factor, and G(U — T) 
would get stronger constraints as a result. Moreover, 
the probe of the cosmic string requires a cone-like and 
degree by degree survey with enormous galaxy sample in 
each degree 2 to increase signal-to- noise ratio (S /N) . This 
kind of survey may be realized by LSST, Dark Energy 
Survey (DES), etc. 

4 Conclusions 

In this paper, we have studied the weak lensing signal 
caused by an infinite cosmic string especially by a wig- 
gly string. In the beginning, we utilize a sensitive indi- 
cator of gravitational distortion, r/c, to replace the non- 
constant distortion angle. We have studied the straight 
jet segment and get the probability to detect a cosmic 
string with current sample and current ability of radio 
telescope. Moreover, the shape of r\Q is different from 
those of other effects and could be distinguished from 
various sources. It is impressive that tjq is constant for 
a straight jet, and r\Q looks like a quasi-sinusoidal curve 
for a spherical lens. This special phenomenon could be 
used to detect a cosmic string. t\q is in the order of 
10 -6 , but A9 and A<j> play an important part to in- 
crease the strength of tjq. Thus, the position of a jet 
determines the value range of rjc] & n d the inclined angle 
of r\Q determines the specific value of t\q. After simple 
calculation, we can determine the relationship between 
distribution of jets number density with rjc through sim- 
ulations. It is the first time that the radio jet is employed 
as an extended gravitational lensing probe to find a wig- 
gly string. It also supplies an independent method to 



study the intrinsic properties of the cosmic string. Due 
to the weak signal of the cosmic string, large scale sur- 
veys is required, and more accurate analysis is needed to 
extract useful signals. It would also be interesting to use 
Expanded Very Large Array (EVLA) to probe the con- 
stant j]g- It would identify the special signal yielded by 
a wiggly string. In addition, we look forward to making 
use of this tool to detect other different kinds of strings. 
As the precision of the radio telescope improved in this 
decade, we would see the appliance of this technique to 
detect a cosmic string and other objects. 

Additionally, ED could also be a sensitive probe of a 
wiggly string. The special area where lensed ED and 
EAD was found would reveal the specific direction of 
a wiggly cosmic string. To be specific, a wiggly string 
would make galaxies more elliptical when they lie in the 
direction of the string. When S = S opt , the SNR of the 
ED reaches the largest value. This optimal mode is ap- 
propriate for us to use the data from SDSS to abstract 
string's signal. For current survey, spirals have higher 
SNR than ellipticals do. With the correction of extinc- 
tion of spirals, we will divide the SDSS sample into small 
groups according to their magnitudes in order to decrease 
the systematic errors in the near future. However, lim- 
ited by the low number density of current sample, EAD 
would not be appropriate for probing a cosmic string. 
In the long run, EAD would be better for cosmic string 
probing because of its independence with ED models. 
Moreover, future data release would make more geomet- 
ric modes of the sky available for putting strong con- 
straints on the abundance of string strongly. 

To sum up, all of these probes require accurate mea- 
surements of a large sample of background objects, and 
ED/EAD requires special geometry of the sky for cur- 
rent data processing. Because ED is easily affected by 
other factors, such as dust extinction and reddening, di- 
viding the galaxy sample into sub-samples is necessary 
for future research on cosmic string detection. It is also 
important to study the weak lensing effect of different 
kinds of cosmic strings and loops. Further, the shear 
and convergence fields of cosmic strings need to be care- 
fully studied too. The next step is to use the sample of 
SDSS to get ED and EAD in each small sky area. With 
more rigorous statistical analyses of various errors, the 
constraints from the real data would be strong enough 
to give upper limit of G(T — U)/c 2 . 



Acknowledgement 

Thanks to Y.F. Huang and Shi Qi for helpful input. 
I am also indebted to H. Brandenberger, P.P. Kronberg 
and the anonymous reviewer for many valuable sugges- 
tions. This work was supported by the National Natural 



X 



Fa-Bo Feng, Radio Jets and Galaxies as Cosmic String Probes, Front. Phys. 




Fig. 6 Left panels show SNRs in the first bin as functions of 6 for different surface number densities and different 
EDs. Middle and right panels show histograms of lensed EDs and EADs with S opt . Each 5 op t is given by the 
left corresponding graph. The error bars in the first and the third (the second and the fourth) histograms are 
generated by choosing the point with its parameters obeying the above distributions about n — 90 A# 2 180 2 /n 2 
(n = 1000 A9 2 /ir) times and repeating this procedure 100 times. The spirals and ellipticals are created through 
100 x n random choices of A9, A<f) e [—0.55, 0.55] and {£, £} according to the ellipticity distributions by Padilla and 
Strauss [38], and randomly choosing the points on an unit sphere which gives the distributions of a, /3,7- 



Fa-Bo Feng, Radio Jets and Galaxies as Cosmic String Probes, Front. Phys. 



xi 



Science Foundation of China (Grant No. 10625313) and gram, Grant 2009CB824800). 
the National Basic Research Program of China(973 Pro- 



References 

[1] T. W. B. Kibble, "Topology of cosmic domains 



and strings," Journal of Physics A Mathematical 



General 9 (Aug., 1976) 1387-1398 



[2] A. Vilenkin and E. P. S. Shellard, Cosmic Strings 
and Other Topological Defects. Cambridge, UK: 
Cambridge University Press, July, 2000. 

[3] J. Polchinski, "Introduction to Cosmic F- and 
D-Strings," ArXiv High Energy Physics - Theory 



e-prints (Dec, 2004) , arXiv:hep-th/0412244 



[4] A. Avgoustidis, "Cosmic String Evolution in Brane 
Inflation," | Journal of Physics C onferenc e Series 



68 no. 1, (May, 2007) 012044-+ 



[5] J. A. Peacock, Cosmological Physics. Cambridge, 
UK: Cambridge University Press, Jan., 1999. 

[6] C. Feng, X. Gao, M. Li, W. Song, and Y. Song, 
"Reheating and cosmic string production," \Nuclear\ 



Physics B 800 (Sept., 2008) 190-203 



arXiv: 0707. 0908 [hep-th] 



[7] B. Chen, Y. Wang, W. Xue, and 

R. Brandenberger, "String Gas Cosmology and 
Non-Gaussianities," ArXiv e-prints (Dec, 2007) , 
|arXiv: 0712 .2477 [hep-th] 



[8] R. Brandenberger, Y. Cai, W. Xue, and X. Zhang, 
"Cosmic Ray Positrons from Cosmic Strings," 



ArXiv e-prints (Jan., 2009) , |arXiv : 0901 . 3474 
| [hep- ph] . 

[9] L. Fang, "Cosmic strings and large scale structure 
of the universe," \Vistas in Astronomy 31 (1988) 
1531-5381 

[10] A. Achucarro and C. J. A. P. Martins, "Cosmic 
Strings," ArXiv e-prints (Nov., 2008) , 
larXiv:0811. 12771 

[11] A. Vilenkin, "Gravitational field of vacuum 



domain walls and strings," Phys. Rev. D 23 (Feb., 
1981) 852-8571 



[12] I. J. R. Gott, "Gravitational lensing effects of 



vacuum strings - Exact solutions," Astrophy s. J. 



288 (Jan., 1985) 422-427 



[13] X. Wu, "Gravitational lenses in the universe.," 
Progress in Astronomy 8 (Mar., 1990) 3-19. 



[14] N. Bevis, M. Hindmarsh, M. Kunz, and 
J. Urrestilla, "Fitting Cosmic Microwave 
Background Data with Cosmic Strings and 
Inflation," Physical Review Letters 100 no. 2, 



|(Jan., 200 8) 02130 1-+[ [arXiv : a stro-ph/0702223 



[15] ct al. B. P. Abbott, "First LIGO search for 
gravitational wave bursts from cosmic 



(super)strings," \Phys. Rev. D 80 no. 6, (Sept., 
|2009) 062002-+I |arXiv : 0904 . 4718| 

[16] L. L. Cowie and E. M. Hu, "The formation of 
families of twin galaxies by string loops," 



Astrophys. J., Lett. 318 (July, 1987) L33-L38 



[17] K. J. Mack, D. H. Wesley, and L. J. King, 

"Observing cosmic string loops with gravitational 



lensing surveys,' 


Phys. Rev. D 76 no. 12, (Dec, 


2007) 123515-+ 


arXiv : astro-ph/0702648 . 



[18] X. Wu, "Lensing effects of a straight cosmic 

string.," Kexue Tongbao 34 (Apr., 1989) 657-660. 

[19] K. Kuijken, X. Siemens, and T. Vachaspati, 



'Microlensing by cosmic strings," \Mon. Not. R. 



Astron. Soc. 384 (Feb., 2008) 161-164 



larXiv: 0707 .29711 

[20] S. Dyda and R. H. Brandenberger, "Cosmic 

Strings and Weak Gravitational Lensing," ArXiv 
e-prints (Oct., 2007) , |arXiv:0710. 1903| 

[21] D. B. Thomas, C. R. Contaldi, and J. Magueijo, 
"Rotation of Galaxies as a Signature of Cosmic 
Strings in Weak Lensing Surveys," Physi cal Review\ 



Letters 103 no. 18, (Oct., 2009) 181301-+ 



arXiv: 0909. 2866 

[22] C. C. Dyer and E. G. Shaver, "On the rotation of 
polarization by a gravitational lens," Astrophys. J., 
\Lett. 390 (May, 1992) L5-L7| 

[23] P. P. Kronberg, C. C. Dyer, E. M. Burbidge, and 
V. T. Junkkarincn, "A technique for using radio 
jets as extended gravitational lensing probes," 
\Astrophys. J., Lett. 367 (Jan., 1991) L1-L4} 

[24] P. P. Kronberg, C. C. Dyer, and H. Roeser, 

"Estimates of the Global Masses of Two Distant 
Galaxies Using a New Type of Astrophysical Mass 



"Laboratory"," Astrophys. J. 472 (Nov., 1996) 



Xll 



Fa-Bo Feng, Radio Jets and Galaxies as Cosmic String Probes, Front. Phys. 



[25] C. R. Burns, C. C. Dyer, P. P. Kronberg, and 

H. Roser, "Theoretical Modeling of Weakly Lensed 



Polarized Radio Sources," Astrophys. J. 613 (Oct. 



2004) 672-681 , arXiv : astro-ph/0406400 



[26] C. R. Burns, Gravitational Lensing of Polarized 
Sources. PhD thesis, University of Toronto, 2002. 

[27] R. A. Pcrlcy, A. H. Bridle, and A. G. Willis, 

"High-resolution VLA observations of the radio jet 
in NGC 6251," \Astrophys. J., Suppl. Ser. 54 (Feb., 



1984) 291-334 



[28] 



N. E. B. Killccn, G. V. Bicknell, and R. D. Ekers, 
"The radio galaxy IC 4296 (PKS 1333 - 33). I. 
Multifrequency Very Large Array observations," 
\Astrophys. J. 302 (Mar., 1986) 306-336 



[29] A. H. Bridle, R. A. Perley, and R. N. Henriksen, 
"Collimation and polarization of the jets in 3C 
219," \Astron. J. 92 (Sept., 1986) 534-545| 

[30] J. Binney, "Testing for triaxiality with kinematic 
data," Mon. Not. R. Astron. Soc. 212 (Feb., 1985) 
767-781. 

[31] B. S. Ryden, "The ellipticity of the disks of spiral 
galaxies," The Astrophysical Journal 601 no. 1, 
(2004) 214.|http: 



//stacks . iop . org/0004-637X/601/i=l/a=214 



[32] J. -P. Uzan and F. Bernardeau, "Cosmic strings 
lens phenomenology: General properties of 



distortion fields," Phys. Rev. D 63 no. 2, (Jan., 
|2001) 023004-+I |arXiv:astro-ph/0004105. 



[33] M. Kamionkowski, "Non-Uniform Cosmological 
Birefringence and Active Galactic Nuclei," ArXiv 
e-prints (Apr., 2010) , |arXiv: 1004.3544| 

[34] V. Faraoni, "The rotation of polarization by 



gravitational waves," New Astronomy 13 (Apr., 



2008) 178-181[ |arXiv:0 709 . 0386 



B. M. Gaensler, R. Beck, and L. Feretti, "The 



origin and evolution of cosmic magnetism," New 



Astronomy Reviews 48 (Dec, 20 04) 1003 -1012 
arXiv : astro-ph/0409100 



[36] A. STARK, "TRIAXIAL MODELS OF BULGE 
OF M31," ASTROPHYSICAL JOURNAL 213 
no. 2, (1977) 368-373. 

[37] K. J. Mack, D. H. Wesley, and L. J. King, 

"Observing cosmic string loops with gravitational 



lensing surveys," Phys. Rev. D 76 no. 12, (Dec, 
|2007) 123515-+[ |arXiv : ast ~ ph/0702648 . 

N. D. Padilla and M. A. Strauss, "The shapes of 
galaxies in the Sloan Digital Sky Survey," Mon. 
\Not. R. Astron. Soc. 388 (Aug., 2008) 1321-1334} 
larXiv: 0802 .08771 

Z. Ivezic, J. A. Tyson, R. Allsman, J. Andrew, 
R. Angel, and for the LSST Collaboration, "LSST: 
from Science Drivers to Reference Design and 
Anticipated Data Products," ArXiv e-prints (May, 
2008) ,|arXiv: 0805 .2366. 



